Dark matter in galaxies: the dark matter particle mass is about 2 keV 
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Warm dark matter (WDM) means DM particles with mass m in the keV scale. For large scales, 
for structures beyond ~ 100 kpc, WDM and CDM yield identical results which agree with observa- 
tions. For intermediate scales, WDM gives the correct abundance of substructures. Inside galaxy 
cores, below ~ 100 pc, iV-body classical physics simulations are incorrect for WDM because at 
such scales quantum WDM effects are important. Quantum calculations (Thomas-Fermi approach) 
provide galaxy cores, galaxy masses, velocity dispersions and density profiles in agreement with the 
observations. All evidences point to a dark matter particle mass around 2 keV. Baryons, which 
represent 16% of DM, are expected to give a correction to pure WDM results. The detection of 
the DM particle depends upon the particle physics model. Sterile neutrinos with keV scale mass 
(the main WDM candidate) can be detected in beta decay for Tritium and Renium and in the 
electron capture in Holmiun. The sterile neutrino decay into X rays can be detected observing DM 
dominated galaxies and through the distortion of the black-body CMB spectrum. The effective 
number of neutrinos, N e ff measured by WMAP9 and Planck satellites is compatible with one or 
two Majorana sterile neutrinos in the eV mass scale. The WDM contribution to N e ff is of the order 
r \ , ~ 0.01 and therefore too small to be measurable by CMB observations. 

. • So far, not a single valid objection arose against WDM. 
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I. INTRODUCTION 

81 % of the matter of the universe is dark. Dark matter (DM) is the dominant component of galaxies. DM interacts 
through gravity. 



* Based on Lectures given by H J de V at NuMass 2013, Milano-Bicocca, February 2013; at Cosmic Frontiers, SLAC, March 2013 and by 

H J de V and N G S at the Chalonge Torino Colloquium 2013, April 2013 
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DM interactions other than gravitational have been so far unobserved, such possible couplings must be very 
weak: much weaker than weak interactions in particle physics. DM is outside the standard model of particle physics. 

The main proposed candidates for DM are: Neutrinos (hot dark matter) back in the 1980's with particle mass m ~ 1 
eV already ruled out, Cold Dark Matter (CDM), weak interacting massive particles (WIMPS) in supersymmetric 
models with R-parity, with particle mass m ~ 10 — 1000 GeV seriously disfavoured by galaxy observations, and finally 
Warm Dark Matter (WDM) , mainly sterile neutrinos with particle mass m ~ 1 keV. 

DM particles decouple due to the universe expansion, their distribution function freezes out at decoupling. The 
characteristic length scale after decoupling is the free streaming scale (or Jeans' scale). Following the DM 
evolution since ultrarelativistic decoupling by solving the linear Boltzmann-Vlasov equations yields (see for example 

0), 

r^on keV /100\* 
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where g ( i equals the number of UR degrees of freedom at decoupling. 

DM particles can freely propagate over distances of the order of the free streaming scale. Therefore, structures at 
scales smaller or of the order of rj eans are erased for a given value of m. 

The observed size of the DM galaxy substructures is in the ~ 1 — 100 kpc scale. Therefore, eq. ljl.ip indicates that 
m should be in the keV scale. That is, Warm Dark Matter particles. This indication is confirmed by the phase-space 
density observations in galaxies pi| and further relevant evidence from galaxy observations [3, 0, ll2T - fl5l ll7Hl9l ] . 

For CDM particles with m ~ 100 GeV we have rj eans ~ 0.1 pc. Hence CDM structures keep forming till scales 
as small as the solar system. This result from the linear regime is confirmed as a robust result by A-body CDM 
simulations. However, it has never been observed in the sky. 

Adding baryons to CDM does not cure this serious problem. There is over abundance of small structures in CDM 
and in CDM+baryons (also called the satellite problem). 

CDM has many serious conflicts with observations as: 



• 



• 



Galaxies naturally grow through merging in CDM models. Observations show that galaxy mergers are rare 
(<10%). 

Pure-disk galaxies (bulgeless) are observed whose formation through CDM is unexplained. 

CDM predicts cusped density profiles: p(r) ~ 1/r for small r. Observations show cored profiles: p{r) bounded 
for small r. Adding by hand strong enough feedback from baryons in the CDM models can eliminate cusps but 
spoils the star formation rate. 

Structures in the Universe as galaxies and cluster of galaxies form out of the small primordial quantum fluctuations 
originated by inflation just after the big-bang. 

The linear small primordial fluctuations grow due to gravitational (Jeans) unstabilities and then classicalize. Struc- 
tures form through non-linear gravitational evolution. Hierarchical structure formation starts from small scales and 
develops to large scales. 

A-body CDM simulations fail to produce the observed structures for small scales less than some kpc. 

Both A-body WDM and CDM simulations yield identical and correct structures for scales larger than some kpc. 

At intermediate scales WDM give the correct abundance of substructures [15| • 

Inside galaxy cores, below ~ 100 pc, A-body classical physics simulations are incorrect for WDM because quantum 
effects are important in WDM at these scales. WDM predicts correct structures for small scales (below kpc) when its 
quantum nature is taken into account |3j. 

The first ingredient in structure formation is the primordial power spectrum P(k). We plot P(k) in fig. [T]for CDM 
and for several examples of WDM. CDM and WDM give identical results for the CMB fluctuations spectrum, this 
spectrum corresponds to large scales > 1 Mpc, in which WDM and CDM coincide. 

All searches of CDM particles (wimps) look for m > 1 GeV |2(j. The fact that the DM mass is in the keV scale 
explains why no detection has been reached so far. Moreover, past, present and future reports of signals of such CDM 




FIG. 1: log 10 P(k) vs. log 10 [fc Mpc h] for CDM in red, for 1 keV WDM in blue, 2 keV in violete and 4 keV WDM in light-blue 
DM particles decoupling in thermal equilibrium. 1 keV WDM sterile neutrinos decoupling out of thermal equilibrium are 
plotted in green. For WDM P(k) is cutted-off on small scales compared with CDM, that is r < 100 (keV/m) 4 / 3 kpc. 



experiments cannot be due to DM detection because the DM particle mass is in the keV scale. The inconclusive 
signals in such experiments should be originated by phenomena of other kinds. Notice that the supposed wimp 
detection signals reported in refs. [2(| contradict each other supporting the notion that these signals are unrelated 
to any DM detection. 

Positron excess in cosmic rays are unrelated to DM physics but to astrophysical sources and astrophysical mecha- 
nisms and can be explained by them [21] . 

Warm dark matter has been the central subject of attention of recent Chalonge Colloquiums whose 'Highlights and 
Conclusions' are online 1171. 



II. QUANTUM PHYSICS IN GALAXIES 



In order to determine whether a physical system has a classical or quantum nature one has to compare the average 
distance between particles with their de Broglie wavelength. 



The de Broglie wavelength of DM particles in a galaxy can be expressed as 

A - h 

*dB — j 



(2.1) 



where v is the velocity dispersion, while the average interparticle distance d can be estimated as 
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where p^ is the average density in the galaxy core. We can measure the classical or quantum character of the system 
by considering the ratio 
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By using the phase-space density, 



and eas. (|2.1[) - (l2.2l) . 1Z can be expressed as Q 
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Notice that 1Z as well as Qh are invariant under the expansion of the universe because the lengths XdB and d both 
scale with the expansion scale factor. 1Z and Qh evolve by nonlinear gravitational relaxation. 
Using now the observed values of Qh from Table U yields 1Z in the range 

2x10-3 f^f <ft<1.4(^f (2.4) 

\ m J \ m J 

The larger value of 1Z is for ultracompact dwarfs while the smaller value of 1Z is for big spirals. 

The ratio 7Z around unity clearly implies a macroscopic quantum object. Notice that 1Z expresses solely in terms 
of Q and hence (H 3 Q/m A ) measures how quantum or classical is the system, here, the galaxy. Therefore, eq. (|2.4p 
clearly shows solely from observations that compact dwarf galaxies are natural macroscopic quantum objects for 
WDM @. 

We see from cq. (f2l| that for CDM, that is for m > GeV, 

Kcdm < 1(T 8 
and therefore quantum effects are negligible in CDM. 

A. WDM Quantum pressure vs. gravitational pressure in compact galaxies 

For an order-of-magnitude estimate, let us consider a halo of mass M and radius R of fermionic matter. Each 
fermion can be considered inside a cell of size Ax ~ l/n3 and therefore has a momentum 

h * i 
p ~ — — ~ a n3 . 

Ax 
The associated quantum pressure P q (flux of the momentum) has the value 

?- - ^ - 
P q = nvp~nvn3 — — n3 . (2.5) 
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where v is the mean velocity given by 



The number density can be estimated as 
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and we obtain from eq. (|2.5[) the quantum pressure 

H 2 ( 3 M 



Pq ~ m R5 {iirm) • (2 - 6) 

On the other hand, as is well known, galaxy formation as all structure formation in the Universe is driven by 
gravitational physics. The system will be in dynamical equilibrium if the quantum pressure is balanced by the 
gravitational pressure 

G M 2 1 

Pa = gravitational force/area = — x (2.7) 

G B ' R 2 Att R 2 K J 

Equating P q = Pq from eas. (|2.7[) - (|2.6l) yields the following expressions for the size R and the velocity v in terms of 
the mass M of the system and the mass m of the particles [3j : 
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Notice that the values of M, R and v are consistent with the observed values of dwarf galaxies. Namely, for M of the 
order 10 6 M Q (which is a typical mass value for dwarf galaxies), R and v give the correct order of magnitude for the 
size and velocity dispersion of dwarf galaxies as displayed in Table U for WDM particle mass in the keV scale. 
These results back the idea that dwarf galaxies are supported by the fermionic WDM quantum pressure eq. 
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TABLE I: Observed values r^, velocity dispersion v, \/Qh, p(0) and Mu covering from ultracompact galaxies to large spiral 
galaxies from refs.[J-[y|. The phase space density is larger for smaller galaxies, both in mass and size. Notice that the phase 
space density is obtained from the stars velocity dispersion which is exp_ected to be smaller than the DM velocity dispersion. 
Therefore, the reported Qh are in fact upper bounds to the true values 



III. QUANTUM FERMIONIC WDM GIVES THE CORRECT GALAXY PROPERTIES AND GALAXY 

PROFILES 



We treat here the self-gravitating fermionic DM in the Thomas-Fermi approximation. In this approach, the central 
quantity to derive is the DM chemical potential /x(r), the chemical potential being the free energy per particle [22j . 
We consider a single DM halo in the late stages of structure formation when DM particles composing it are non- 
relativistic and their phase-space distribution function /(i,r,p) is relaxing to a time-independent form, at least for 
r not too far from the halo center. In the Thomas-Fermi approach such a time-independent form is taken to be a 
energy distribution function f(E) of the conserved single-particle energy E — p 2 /(2m) — /i, where m is the mass of 
the DM particle and \x is the chemical potential 



M(r) = /Uo-m^(r) 

with <p(v) the gravitational potential and fj,o some constant. We consider the spherical symmetric case. 
Here, the Poisson equation for <j)(r) is a nonlinear and selfconsistent equation 
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where the mass density p(r) is a function of p(r) and G is Newton's constant. p(r) is expressed here as a function 
of p(r) through the standard integral of the DM phase-space distribution function over the momentum for Dirac 
fermions as 

pM = ^fW /(£-<*•>), (3.3) 

Another standard integral of the DM phase-space distribution function is the pressure 

p «-3^r dp// (£- M(r) )- (3 - 4) 

From p(r) and P(r) other quantities of interest, such as the velocity dispersion a(r) and the phase-space density Q(r) 
can be determined as 

° 2 (r) = ^ , Q(r) = 4±. (3.5) 

p(r) G a (r) 

We see that p{r) fully characterizes the fermionic DM halo in this Thomas-Fermi framework. The chemical potential 
is monotonically decreasing in r since eq. (I3.2[) implies 
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Af(r)=4vr/ dr' r' 2 p(r') . (3.6) 
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Moreover, the fermionic DM mass density p is bounded at the origin due to the Pauli principle J3J, and therefore the 
proper boundary condition at the origin is 

^(0)=0. (3.7) 

Eqs. (|3.2p and (|3.3p provide an ordinary nonlinear differential equation that determines selfconsistently the chemical po- 
tential p,(r) and constitutes the Thomas-Fermi semi-classical approach. We obtain a family of solutions parametrized 
by the value of po = A*(0) [3]. 

In this semi-classical framework the stationary energy distribution function f(E) must be assigned beforehand. In 
a full-fledged treatment one would solve the cosmological WDM evolution since decoupling till today, including the 
quantum WDM dynamics in the evolution which become important in the non-linear stage and close enough to the 
origin. 

We integrate the Thomas-Fermi nonlinear differential equations (|3.2[) - (|3.3[) from r = till the boundary r = R = 
i?200 ~ Rvir defined as the radius where the mass density equals 200 times the mean DM density [3j . 
We define the core size r^ of the halo by analogy with the Burkert density profile as 

- i — ^ = t , r h = Z & • (3.8) 

Po 4 

where po = p(0) and Iq is the characteristic length that emerges from the dynamical equations 
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To explicitly solve eqs. (|3.2[ )- (l3.3p we need to specify the distribution function ^{E/ Eq). But many important proper- 
ties of the Thomas-Fermi semi-classical approximation do not depend on the detailed form of the distribution function 
^(E/Eq). Indeed, a generic feature of a physically sensible one-parameter form ^{E/Eq) is that it should describe 
degenerate fermions for Eq — > 0. That is, ^(E/Eq) should behave as the step function 9(—E) in such limit. In the 
opposite limit, p/ Eq — > — oo, ^>(E/Eq) describes classical particles, namely a Boltzmann distribution. As an example 
of distribution function, we consider the Fermi-Dirac distribution 

**°WEo) = ^rr x ■ (3-10) 

We define the dimensionless chemical potential i/(r) as 

v(r) = p{t)/Eq and vq = p(0)/E o . 



Large positive values of the chemical potential at the origin i/g 3> 1 correspond to the degenerate fermions limit which 
is the extreme quantum case, and oppositely, Vq <C — 1 gives the diluted regime which is the classical limit. In this 
classical regime the Thomas-Fermi equations (|3.2|) - (|3.3|) become exactly the equations for a self-gravitating Boltzmann 
gas. 

We display in fig. [2]the density and velocity profiles Q. Namely, we plot p(r)/po and cr(r)/cr(0) as functions oir/R 
for Vq = z/(0) = —5, 0, 5, 15, 25 and for the degenerate fermion limit Vq — > +oo. The obtained fermion profiles are 
always cored. The sizes of the cores r h defined by eq. (|3.8[) are in agreement with the observations, from the compact 
galaxies where r^ ~ 35 pc till the spiral and elliptical galaxies where r^ ~ 0.2 — 60 kpc. The larger and positive is vq, 
the smaller is the core. The minimal core size arises in the degenerate case v n — > +oo (compact dwarf galaxies). 
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FIG. 2: Density and velocity profiles, p(r)/po and cr(r)/er(0), as functions of r/lo for different values of the chemical potential 
at the origin vq [3(]. Large positive values of vo correspond to compact galaxies, negative values of vq correspond to the classical 
regime describing spiral and elliptical galaxies. All density profiles are cored. The sizes of the cores r^ defined by eg. 1)3.81) are 
in agreement with the observations, from the compact galaxies where rn ~ 35 pc till the spiral and elliptical galaxies where 
ru ~ .2 — 60 kpc. The larger and positive is vo, the smaller is the core. The minimal one arises in the degenerate case uq — > +cxj 
(compact dwarf galaxies). 



In the left panel of fig. [3] we plot the dimensionless quantity Q 

(keW O(0) • 
In the right panel of fig. 02 we plot instead the dimensionless product 



M h 1 M 
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where Mh is the halo mass, namely the galaxy mass inside the core radius r/, defined by eq. 
consider the two values m — 1 and 2 keV and we put in the abscissa the product 
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FIG. 3: In the left panel we display the galaxy phase-space density ft 3 Q(0)/(keV) 4 obtained from the numerical resolution 



of the Thomas-Fermi eqs. (I3.2I) - (|3.3[) for WDM fermions of mass m = 1 and 2 keV versus the ordinary logarithm of the 
product log 10 {rh [pc 3 po/Mq]®} in parsecs Q|. The red stars * are the observed values of ft 3 Q(0)/(keV) 4 from Table[I] Notice 
that the observed values Qh from the stars' velocity dispersion are in fact upper bounds for the DM Qh and therefore the 
theoretical curve is slightly below them. In the right panel we display the galaxy mass (M /Mq)-^ 'M@/[po pc 3 ] obtained from 
the numerical resolution of the Thomas- Fermi eqs. (|3.2|) - (|3.3l) for WDM fermions of mass m — 1 and 2 keV versus the product 

Th. [pc 3 po/Mq]® in parsecs |g]. The red stars * are the observed values of (M/Mq)\JMq/\pq pc 3 ] from Table [1] Notice that 
the error bars of the observational data are not reported here but they are at least about 10 — 20%. 



where r^ is the core radius. The phase-space density Q(0) and the galaxy mass M^ are obtained by solving the Thomas- 
We have also superimposed the observed values H 3 Q/j/(keV) 4 and M/ja/Mq/^o pc 3 ] (m/keV) 



Fermi eqs. (|32D-(|33D 

from Table UJ Notice that the observed values Qh from the stars' velocity dispersion are in fact upper bounds for the 
DM Qh. This may explain why the theoretical Thomas-Fermi curves in the left panel of fig. [3] appear slightly below 
the observational data. Notice also that the error bars of the observational data are not reported here but they are 
at least about 10 - 20%. 

The phase space density decreases from its maximum value for the compact dwarf galaxies corresponding to the 
limit of degenerate fermions till its smallest value for large galaxies, spirals and ellipticals, corresponding to the 
classical dilute regime. On the contrary, the halo radius r/j and the halo mass Mh monotonically increase from the 
quantum (small and compact galaxies) to the classical regime (large and dilute galaxies). 

Thus, the whole range of values of the chemical potential at the origin v from the extreme quantum (degenerate) 
limit i/o > 1 to the classical (Boltzmann) dilute regime vq <C — 1 yield all masses, sizes, phase space densities 
and velocities of galaxies from the ultra compact dwarfs till the larger spirals and elliptical in agreement with the 
observations (see Table fl}. 

From figs. [3] we can extract important information on the fermion particle WDM mass. 

We see from the left panel fig. [3]that decreasing the DM particle mass m moves the theoretical curves ft 3 Q(0)/(kcV) 4 
towards smaller Q(0) values and larger galaxy sizes, one over each other. In the right panel of fig. [3] we see that 
decreasing the DM particle mass m displaces the theoretical curves (Mh/Mo)^J Mq/[p q pc 3 ] towards larger galaxy 
masses and sizes, one over each other. 



The small galaxy endpoint of the curves in figs. [3] corresponds to the degenerate fermion limit vq — > +00 and its 
value depends on the WDM particle value to. For increasing to, the small galaxy endpoint moves towards smaller 
sizes while for decreasing to, it moves towards larger sizes. 

We see from figs. [3] that decreasing the particle mass beyond a given value, namely for particle masses to < 1 
keV, the theoretical curves do not reach the more compact galaxy data. Therefore, to < 1 keV is ruled out as WDM 
particle mass. 

For growing to > keV the left part of the theoretical curves corresponding to the lower galaxy masses and sizes, will 
not have observed galaxy counterpart. Namely, increasing to 3> keV would show an overabundance of small galaxies 
(small scale structures) which do not have counterpart in the data. This is a further indication that the WDM particle 
mass is approximately around 2 keV in agreement with earlier estimations [7J, llCf . 

The overabundance of small scale structures that appears here for to ^> keV is a consequence only of the DM 
particle mass value. This result is independent of the cosmic evolution, namely the primordial power spectrum and 
the structure formation dynamics. In particular, this argument rules out hybrid Wimp models including those with a 
cutoff in the primordial power spectrum for small scales (as in WDM) but with a DM particle mass in the GeV scale 
or above [11] . 

In addition, the galaxy velocity dispersions turn to be fully consistent with the galaxy observations in Table [J Q . 

Adding baryons to CDM simulations have been often invoked to solve the serious CDM problems at small scales. It 
must be noticed however that the excess of substructures in CDM happens in DM dominated halos where baryons are 
especially subdominat and hence the effects of baryons cannot drastically modify the overabundance of substructures 
of the pure CDM results. 

The influence of baryon feedback into CDM cusps of density profiles depends on the strength of the feedback. For 
normal values of the feedback, baryons produce adiabatic contraction and the cusps in the density profiles become 
even more cuspy. 

Using the baryon feedback as a free parameter, it is possible to exagerate the feedback such to destroy the CDM 
cusps but then, the star formation ratio disagrees with the available and precise astronomical observations. Moreover, 
" semi- analytic (CDM + baryon) models" have been introduced which are just empirical fits and prescriptions to some 
galaxy observations. 

In addition, there are serious evolution problems in CDM galaxies: for instance pure-disk galaxies (bulgeless) are 
observed whose formation through CDM is unexplained. 

In summary, adding baryons to CDM simulations bring even more serious discrepancies with the set of astronomical 
observations. 

In this paper spherical symmetry is considered for simplicity to determine the essential physical galaxy properties 
as the classical or quantum nature of galaxies, compact or dilute galaxies, the phase space density values, the cored 
nature of the mass density profiles, the galaxy masses and sizes. It is clear that DM halos are not perfectly spherical 
but describing them as spherically symmetric is a first approximation to which other effects can be added. In ref. 
[9] we estimated the angular momentum effect and this yields small corrections. The quantum or classical galaxy 
nature, the cusped or cored nature of the density profiles in the central halo regions can be captured in the spherically 
symmetric treatment. 

Our spherically symmetric treatment captures the essential features of the gravitational dynamics and agree with 
the observations. Notice that we are treating the DM particles quantum mechanically through the Thomas-Fermi 
approach, so that expectation values are independent of the angles (spherical symmetry) but the particles move and 
fluctuate in all directions. Namely, this is more than treating purely classical orbits for particles in which only radial 
motion is present. The Thomas-Fermi approach can be generalized to describe non-spherically symmetric and non- 
isotropic situations, by considering distribution functions which include other particle parameters like the angular 
momentum. 

To conclude, the galaxy magnitudes: halo radius, galaxy masses and velocity dispersion obtained from the Thomas- 
Fermi quantum treatment for WDM fermion masses in the keV scale are fully consistent with all the observations for 
all types of galaxies (see Table [J). Namely, fermionic WDM treated quantum mechanically, as it must be, is able to 
reproduce the observed DM cores and their sizes in galaxies Q. These results strenght the discussion in sec. Ill Al that 
compact galaxies are supported against gravity by the fermionic WDM quantum pressure. 

It is highly remarkably that in the context of fermionic WDM, the simple stationary quantum description provided 
by the Thomas-Fermi approach is able to reproduce such broad variety of galaxies. 

Baryons have not yet included in the present study. This is fully justified for dwarf compact galaxies which are 
composed today 99.99% of DM. In large galaxies the baryon fraction can reach values up to 1 - 3 %. Fermionic WDM 
by itself produces galaxies and structures in agreement with observations for all types of galaxies, masses and sizes. 
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Therefore, the effect of including baryons is expected to be a correction to these pure WDM results, consistent with 
the fact that dark matter is in average six times more abundant than baryons. 

IV. WDM GIVES THE CORRECT ABUNDANCE OF SUBSTRUCTURES 

It is known since some time through A-body simulations that WDM alleviates the CDM satellite problem [12|, LL3[ 
and the CDM voids problem fl~H ]. 

WDM subhalos turns to be less concentrated than CDM subhalos. WDM subhalos have the right concentration to 
host the bright Milky Way satellites [l3J]. 

The ALFALFA survey has measured the velocity widths in galaxies from the 21cm HI line. This is a test for 
substructure formation. The confrontation of the ALFALFA survey with the substructures from A-body simulations 
clearly favours WDM over CDM 15]. A particle mass around ~ 2 keV is favoured by the ALFALFA survey. 

In summary, WDM produces the correct substructure abundance at zero redshift. 

Data on galaxy substructure for redshift z < 10 becomes now available. 

In ref. [18|] the evolution of the observed AGN luminosity function for 3 < z < 6 is contrasted with the WDM and 
CDM simulations. WDM is clearly favoured over CDM by the observational data. 

In ref. [19[ the number of observed structures vs. the star formation rate for z = 5, 6, 7 andf 8 is contrasted with 
the results from WDM and CDM simulations. Again, WDM turns to be clearly favoured by the observations over 
CDM. 

At intermediate scales where WDM and CDM give non-identical results and quantum effects are negligible, iV-body 
classical simulations are reliable. Contrasting such A-body classical simulations results with astronomical observations 
at zero and non-zero redshifts clearly favours WDM over CDM. 

For larger scales > 100 kpc, CDM and WDM A-body classical simulations are reliable and give identical results in 
good agreement with astronomical and cosmological observations. 

V. DETECTION OF KEV MASS STERILE NEUTRINOS 

Sterile neutrinos v s are mainly formed by right-handed neutrinos vr plus a small amount of left-handed neutrinos 
vi,. Conversely, active neutrinos v e are formed by vl plus a small amount of vr- 

v s — VB. + dv L , v a = v L + 9 v R ■ 

The name sterile neutrino were was introduced by Bruno Pontecorvo in 1968. They are singlets under all symmetries 
of the Standard Model of particle physics. Sterile neutrinos do not interact through weak, electro-magnetic or strong 
interactions. 

WDM v s are typically produced in the early universe from active neutrinos through mixing, namely, through a 
bilinear term 9 v s v a in the Lagrangian. 

The appropriate value of the mixing angle 9 to produce enough sterile neutrinos v s accounting for the observed 
total DM depends on the particle physics model and is typically very small: 

9 ~ 10~ 3 - 10~ 4 . 

The smallness of 9 makes sterile neutrinos difficult to detect in experiments. 

Sterile neutrinos can be detected in beta decay and in electron capture (EC) processes when a v s with mass in the 
keV scale is produced instead of an active v a : 

3 H 1 => 3 He 2 + e~ + v e , 187 Re =>■ 187 Os + e~ + v e . 

In beta decays when a v s is produced instead of a v & in the decay products the electron spectrum is slightly modified 
at energies around the v s mass (~ keV). Such event can be inferred observing the electron energy spectrum. A 'kink' 
should then appear around the energy of the v s mass. 



11 

In electron capture processes like: 

163 Ho + e~ =*> 163 Dy* + v e , 

when a sterile neutrino v a with mass in the keV scale is produced instead of an active v e , the observed nonradiative 
de-excitation of the excited Dysprosium Dy* is different to the case where an active v e shows up. 

The available energies for these beta decays and EC are 

Q( 187 Re) = 2.47 keV , Q( 3 #i) = 18.6 keV , Q( 163 Ho) ~ 2.5 keV. (5.1) 

In order to produce a sterile neutrino with mass m, Q must be larger than to. However, in order to distinguish the 
sterile neutrino v s from a practically massless active neutrino v a , Q must be as small as possible. This motivates the 
choice of the nuclei with the lowest known Q in eq. (|5.1|) . 

For a theoretical analysis of v s detection in Rhenium and Tritium beta decay see ref. [23j and references therein. 
Present experiments searching the small active neutrino mass also look for sterile neutrinos in the keV scale: 

• MARE (Milan, Italy), Rhenium 187 beta decay and Holmiun 163 electron capture [28| . 

• KATRIN (Karlsruhe, Germany), Tritium beta decay [I<|. 

• ECHo (Heidelberg, Germany), Holmiun 163 EC [13]. 

• Project 8 (Seattle, USA), Tritium beta decay J3lj . 
The more popular sterile neutrino models nowadays are: 

• The Dodelson-Widrow (DM) model (1994): sterile neutrinos are produced by non-resonant mixing from active 
neutrinos. 

• The Shi- Fuller model (SF) (1998): sterile neutrinos are produced by resonant mixing from active neutrinos. 

• z/MSM model (2005): sterile neutrinos are produced by a Yukawa coupling from the decay of a heavy real scalar 
field x- 

• Models based on: Froggatt-Nielsen mechanism, flavor symmetries, Qe, split see-saw, extended see-saw, inverse 
see-saw, loop mass. Furthermore: scotogenic, LR symmetric, etc. See for a recent review j24j |. 

WDM particles decoupling ultrarelativistically in the first three WDM particle models behave just as if their 
masses were different Q. The masses of WDM particles in the first three models which give the same primordial 
power spectrum can be related according to the formula [2j (FD = thermal fermions): 

— — — ~ 2.85 I -— — J , m S F ^ 2.55 m FD , to^msm ^ 1-9 m FD . 

For the primordial power spectrum, there is a degeneracy between these three models. 

Sterile neutrinos v s decay into active neutrinos v a plus X-rays with an energy to/2 [25j . The lifetime of v s is about 
~ 10 lx x age of the universe. The value of the lifetime depends on the particle physics neutrino model. 

These X-rays may be seen in the sky looking to galaxies [26J. See [27| for a recent review. 
Some future observations of X-rays from galaxy halos: 

• DM bridge between M81 and M82 ~ 50 kpc. Overlap of DM halos. Satellite projects: Xenia (NASA) [Hj]. 
ASTRO-H (Japan) J2J|. 

• CMB: WDM decay distorts the blackbody CMB spectrum. The projected PIXIE satellite mission can measure 
WDM sterile neutrino mass by measuring this distortion [34|. 

Active neutrinos are very abundant in supernovae explosions and in these explosions sterile neutrinos are produced 
too. Hence, bounds on the presence of sterile neutrinos can be obtained contrasting to supernovae observations. The 
results from supernovae do not constrain 9 provided 1 < to < 10 keV [35J. 
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VI. STERILE NEUTRINOS AND CMB FLUCTUATIONS 



CMB fluctuations data provide the effective number of neutrinos, N e ff. This effective number N e ff is related in 
a subtle way to the real number of active neutrinos (three) plus the number of sterile neutrinos with masses much 
smaller than the electron mass m e [36| • 

WDM should decouple early at temperatures beyond the Fermi scale because WDM is not in the Standard Model of 
Particle Physics (SM). WDM couple much weakly than weak interactions to the SM particles. Therefore, the number 
of ultrarelativistic degrees of freedom at decoupling g(Td) includes all SM particles and probably beyond. We have 
9SM = 427/4 for the SM and gmssM = 915/4 for the minimal supersymmetric SM. 

Entropy conservation determines the WDM contribution to the effective number of neutrinos N e s [371 • One WDM 
sterile neutrino decoupling at Td contributes to N e ff at recombination (re) by 



AN 



WDM 



T 

-L rr. 



9r 



-,4/3 



9(T d ) 



(6.1) 



At recombination z = 1090, we have g rc = 29/4 and eq. ffTTj) gives for the SM and the MSSM: 



A atWDM 



0.02771 



AiV 



WDM 

MSSM 



0.01003. 



Such WDM contributions are too small to be measurable by present CMB anisotropy observations. Hence, Planck 
and WMAP results cannot provide information about these sterile neutrino WDM contributions. 
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In addition, Planck results [38| are compatible with one or two Majorana sterile neutrinos in the eV mass scale 
[36} . The possibility of a eV sterile neutrino is important for the whole subject of neutrinos and hence for WDM. 
The existence of one sterile neutrino in the eV scale suggests the existence of further sterile neutrinos with different 
masses, including a keV mass WDM sterile neutrino. 

VII. FUTURE PERSPECTIVES AND STERILE NEUTRINO DETECTION 

WDM particle models must explain the baryon asymmetry of the universe. This is a strong constraint on sterile 
neutrino models which must be worked out for each model. 

Combining particle, cosmological and galaxy results for sterile neutrinos at different mass scales [31 KLOl I39l44l| ■ an 
appealing mass neutrino hierarchy appears: 

• Active neutrino: ~ mili eV 

• Light sterile neutrino: ~ eV 

• Dark Matter sterile neutrino: ~ keV 

• Unstable sterile neutrino: ~ MeV.... 

This scheme may represent the future extension of the standard model of particle physics. 

In order to falsify WDM, comprehensive theoretical calculations showing substructures, galaxy formation and 
evolution including the quantum WDM effects in the dynamical evolution are needed to contrast with the astronomical 
observations. 

In such WDM theoretical calculations the quantum pressure must be necesarily included. These calculations should 
be performed matching the semiclassical Hartree-Fock (Thomas- Fermi) dynamics where the dimensionless phase-space 
density in high enough, namely, ft? Q/m? > 0.1 with the classical evolution dynamics where ft? Q/m? <C 1. These 
are certainly not easy numerical calculations but they are unavoidable! 

Richard P. Feynman foresaw the necessity to include quantum physics in simulations in 1981 [42] 

"I'm not happy with all the analyses that go with just the classical theory, because nature isn't 

classical, dammit, and if you want to make a simulation of nature, you'd better make it quantum 

mechanical, and by golly it's a wonderful problem, because it doesn't look so easy." 

Sterile neutrino detection depends upon the particle physics model. There are sterile neutrino models where the 
keV sterile is stable and thus hard to detect (see for example [24|V 

Detection may proceed through astronomical observations of X-ray keV sterile decay from galaxy halos and by 
direct detection of steriles in laboratory experiments. 

Mare [28|], Katrin [29[, ECHo J3(| and Project 8 (31J are expected to provide bounds on the mixing angles. However, 
for a particle detection, a dedicated beta decay experiment and/or electron capture experiment seem necessary to 
find sterile neutrinos with mass around 2 keV. In this respect, calorimetric techniques seem well suited. 

The best nuclei for study are: 187 Re and Tritium for beta decay and 163 Ho for electron capture 

The search of DM particles with mass around 2 keV is a promisory avenue for future trascendental discoveries. 
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